Quantum mechanical study of a generic quadratically coupled optomechanical system 
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Typical optomechanical systems involving optical cavities and mechanical oscillators rely on a 
coupling that varies linearly with the oscillator displacement. However, recently a coupling varying 
instead as the square of the mechanical displacement has been realized, presenting new possibilities 
for non-demolition measurements and mechanical squeezing. In this article we present a quantum 
mechanical study of a generic quadratic-coupling optomechanical Hamiltonian. First, neglecting 
dissipation, we provide analytical results for the dressed states, spectrum, phonon statistics and 
entanglement. Subsequently, accounting for dissipation, we supply a numerical treatment using a 
master equation approach. We expect our results to be of use to optomechanical spectroscopy, state 
transfer, wavefunction engineering, and entanglement generation. 

PACS numbers: 42.50.Pq, 42.60.Da, 42.50.Ct,42.50.Dv 



I. INTRODUCTION 

Optomechanics, which involves the interaction be- 
tween high-finesse optical resonators and high quality 
mechanical oscillators, is currently of great interest as a 
new frontier in quantum mechanics [l|-|3|. Developments 
that justify such interest include the measurement of me- 
chanical displacements below the standard quantum limit 
d, Q , the observation of back-action effects Q , the op- 
tical cooling of an oscillator to its quantum mechanical 
ground state [7| and observation of its zero-point mo- 
tion [8|. These pioneering experiments, as well as many 
others, have relied on an optomechanical coupling that 
varies linearly with the displacement of the mechanical 
oscillator. 

Recently however, an optomechanical coupling that 
varies quadratically with the mechanical displacement 
has been experimentally demonstrated in a number of 
configurations, including membranes [9|, |lO[ and ultra- 
cold atoms [ll| in optical cavities, and microdisk res- 
onators [l2l- As opposed to the linearly-coupled sys- 
tems, the quadratic coupling offers new possibilities for 
non-demolition measurements [9|, cooling and squeezing 
[l3l - [la |. optical springs [131, ^-^^ photon transport flSj . 

In this article we present a quantum mechanical study 
of a generic quadratically coupled optomechanical Hamil- 
tonian that can model any of the three experimental real- 
izations mentioned above. Our work provides the dressed 
states of the system, and clarifies further advantages of 
quadratic over linear coupling, such as a non-vanishing 
optomechanical entanglement (in the absence of dissipa- 
tion) , a higher nonlinearity of the spectrum with respect 
to photon number, and the possibility of achieving me- 
chanical squeezing without a projective measurement of 
the system wavefunction [19l - l2l| . We note that some 
of the work presented here is complementary to that of 
Ref. [l7|, which includes a description of oscillator oc- 



cupation probabilities, quadrature squeezing, and back- 
action effects. 

The remainder of this article is arranged as follows. 
Section |lT] considers a quadratic optomechanical system 
in the absence of noise and dissipation, and presents the 
eigenstates, eigenvalues, entanglement, phononic prop- 
erties and nonclassical state engineering. Section HIT] de- 
scribes the inclusion of noise and dissipation using a mas- 
ter equation approach. Section HVl supplies a conclusion. 



II. QUADRATICALLY COUPLED 

OPTOMECHANICAL SYSTEM WITHOUT 

NOISE AND DISSIPATION 

The generic Hamiltonian for a quadratically coupled 
optomechanical system may be written as 



H — Hs + Ho + H„ 



(1) 



where Hs describes the closed system consisting of an 
optical mode and a mechanical resonator, while Ho and 
H„i describe the coupling of each of these modes to the 
external environment. The system Hamiltonian is given 
in detail by [H, (H 



Hs 
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where a{a^) and 6(6') are the annihilation(creation) op- 
erators obeying the bosonic commutation rules [a^a^ = 
[6,6^] = 1, and Wo and uJm are the resonant frequen- 
cies of the optical and mechanical modes respectively. 
The strength of the optomechanical coupling in Eq. ^ 
is measured by g, whose explicit form depends on the 
particular experimental realization [9|, I 111 . Il2l | . An exter- 
nal field, if present, drives the system at frequency w^ 



and is represented by the coupling 

1/2 



E 



hhJn 



(3) 



where 70 is the rate of optical damping and Pin is the 
input power into the optical resonator. 

The derivation of Eq. ^ requires several assumptions, 
as described in detail in Ref. [23] ■ There are two essen- 
tial approximations. First, it is assumed that photons are 
not scattered into other modes of the optical resonator 
by the mechanical oscillator, i.e. Wm is much smaller 
than the free spectral range of the optical cavity. Sec- 
ond, the amplitude of mechanical motion is assumed to 
be much smaller than an optical wavelength, which im- 
plies that only the lowest-order dependence of the res- 
onator frequency on oscillator displacement needs to be 
considered. 

In this section, we neglect environmental effects and 
set Ho = Hm = 0, based on three motivations. First, 
the physics of the closed quantum system modeled by 
the Hamiltonian Hs is already quite rich and interesting. 
Second, it is easier to understand the system behavior in 



the absence of external effects. Third, experiment [24 1 
and theory [25| are both moving towards optomechanical 
configurations in which environmental effects are greatly 
reduced. If optical loss is negligible {Ho = 0), it follows 
that 7o = 0. Also, in the absence of driving, E = Q. In 
that case, Eq. ^ yields [i3| 



Hq = hojQo'a + hujjnb b + hga*a{h* + 6) 



(4) 



Clearly the operator Na = a) a commutes with the Hamil- 
tonian of Eq. Q, and thus the number of photons is 
conserved. In contrast, [6^6, Hq] 7^ 0, showing that 
the number of phonons is not conserved. We note that 
zero-drive Hamiltonians have been used by multiple au- 
thors to analyze optomechanical systems quantitatively 

fiSIIKiilia. 



A. Eigenstates 

In this section we present the dressed states of the sys- 
tem, the knowledge of which is important to applications 
such as state transfer and wavefunction engineering. We 
assume that each eigenstate of the Hamiltonian of Eq. ([4]) 
can be written simply as a product of the optical state 
and the mechanical state. Further, since [a) a^ Hq] = 0, 
we expect the optical contribution to the system eigen- 
state to be simply a number state \n)c, where the sub- 
script c stands for 'cavity'. For g = the mechanical 
contribution is also a number state |fc)m- Although this 
is no longer true if 5 7^ 0, the form of the third term 
in Eq. (|4]) suggests that the optical field 'squeezes' the 



g/w,„ = 0.003 



200 400 600 800 



FIG. 1: The n— dependence of the squeezing parameter r{n) 
of Eq. (0 for g/co^ = 0.003 [H. 



position of the mechanical oscillator in proportion to the 
number of optical quanta n. Hence we may expect a 
'squeezed' number state as a contribution from the me- 
chanical mode. We thus make the ansatz 



\k,n) = S[r{n)]\k)m\n)c, 
for the system eigenstate, where 

b[r(n)\ = e 2 V J ^ 



(5) 



(6) 



is the squeezing operator of quantum optics and r{n) 
is the cavity photon number dependent dimensionless 
squeezing parameter. 

We find that \k,n) of Eq. ([5]) satisfies the eigenvalue 
equation H\k,n) = Ek^n\k,n) if r{n) is a positive real 
number given by 



r{n) = -tanh 



1 + 



2gn 



(7) 



which is plotted as a function of the cavity photon num- 
ber in Fig. [T] for experimentally relevant parameters [13 |. 
As in Fig. [1] since we are primarily interested in the 
behavior of the mechanical oscillator, we scale all fre- 
quencies and times in our problem appropriately by the 
mechanical frequency Wm in the remainder of the article. 

With regard to the eigenstates, we make three obser- 
vations. First, we find that for g = we regain from 
Eq. ([5]) the expected product |n)c|fc)„ for the system 
eigenstate. Second, we note that the dressed states of 
Eq. ^ arc orthonormal, i.e. {k',n'\k,n) = Sk'^k^ri 
unlike the eigenstates provided in previous studies [Tj 
Lastly, we warn the reader that k does not denote the 
number of phonons in an eigenstate, unless n = 0. 

The eigenstates \k,n) are therefore products of optical 
number states and squeezed mechanical number states, 



where the squeezing parameter for each mechanical state 
is a function of the number of photons in the corre- 
sponding optical number state. They are analogous to 
the linear optomechanical coupling eigenstates which are 
products of optical number states and displaced mechan- 
ical number states, where the displacement parameter 
for each mechanical state is a function of the number 
of p hotons in the corresponding optical number state 

Since the properties of number states are very well 
known in quantum optics, we will not describe them here. 
Squeezed number states have also been studied earlier in 
the literature |28|, and we will only make a comment on 
the phonon statistics of the eigenstate. The expectation 
value of the phonon number is given by 

(fc, n\b''b\k, n) ^ k + {2k + 1) sinh^ r{n), (8) 

which depends on both k as well as n, as expected. We 
also find, 

(fc,n| (b^b)' \k,n} = {b^bf + (k' + k + l)^sinh'2rin)^ 

(9) 
Using the above results, we can find analytically the 
phononic Mandel parameter defined by 



_ ((bt&)2) _ (btj,)2 



(10) 



This parameter indicates the nature of the phonon statis- 
tics. It is straightforward to verify that an eigenstate 
|fc,n) of Eq. ([5]) can display any type of phonon statis- 
tics, depending on the parameters. For example, for 
k ^ 0,n = (phononic number state) we find Qk = — 1, 
implying sub-Poissonian statistics. Also, for k = 0,n ^ 
(squeezed mechanical vacuum) it follows that Qk = 
cosh2r(7T,) > 1, which corresponds to super- Poissonian 
phonon statistics. Lastly, it can easily be shown that 
Qk ~ 0, i.e. the eigenstate phonon statistics are Poisso- 
nian, if the eigenstate parameters satisfy 



sinh r(n) 
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B. Eigenvalues 



(11) 



We now discuss the eigenvalues Ek,n of the Hamilto- 
nian of Eq. (jj), which are important to applications such 
as optomechanical spectroscopy [21| and photon blockade 
|26l |27| and are given by 
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FIG. 2: The spectrum of the quadratically coupled optome- 
chanical Hamiltonian without driving, see Eq. (|12|l . For view- 
ing clarity ctJo = '2idm has been used. 
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FIG. 3: The numerically calculated spectrum of the Hamil- 
tonian of Eq. ([2]) with the same parameters as Fig. [2] and 
including driving (E/ujm = 1). It can be seen that some of 
the crossings of Fig. [2] namely those between states differing 
in photon number by one, have turned into avoided crossings 
l/due to the presence of external driving. 



This expression has been derived earlier by Rai et. al 
[l7j . As can be seen in Fig. [21 the spectrum is relatively 
simple for one-photon excitations, but quite rich when 
few-photon excitations are considered. It is interesting 
to ask how the spectrum is changed by the presence of 
driving. Therefore, for comparison, we have also shown a 
numerical calculation of the Floquet spectrum of Eq. ([2]) 
in Fig. [3l which reveals the effects of driving. 

It is worth remarking that the spectrum of the quadrat- 
ically coupled Hamiltonian is anharmonic to all orders in 
n, unlike the case of linear coupling, where the anhar- 
monicity is quadratic in n |2ll. l26l |27| . In fact, the an- 
harmonicity in the present quadratic case is of the same 



functional form as in the well-known Jaynes-Cummings 
spectrum, i.e. of the type {A + Bn)^/^, where A and B 
are quantities independent of the photon number n. The 
spectrum described in this section will be observable in 
the limit of strong single-photon coupling, which is being 
approached in optomechanical systems [26|. 



C. Time evolution of the state vector 

Using the Hamiltonian of Eq. (JU and standard dis- 
entangling techniques from quantum optics [29[, we can 
evaluate the time evolution operator in a frame rotating 
at the optical frequency wg 

jj ^ g«a;„t/2g|(rfc'-«6t')eJ')(6tb+|)^ (13) 



where 



with 



and 



where 



^ = e<"+t)sinh-^('i:^^ 
V X 



r = 2.ga^a, x = [^rn (w,„ + 4ga^a)] 



t„Mi/2 



77 = tan 



^ fptanxt 
2x 



p = —2 [LJm + 2ga^aj 



(14) 



(15) 



(16) 



(17) 



We note that rj is purely real. 

With the time evolution operator of Eq. (fT3|) we 
can, given an initial state vector |'0(O)), find the vector 
\ip{t)) = U\ip{0)) at any later time t. We will assume as 
our initial state 

|V(0)) = |0)„,|a)e, (18) 

where |0)„i is the mechanical ground state, and \a)c is 
an optical coherent state. Such an initial state can be 
realized in experiments which have placed mechanical os- 
cillators in thermal states with occupation number prac- 
tically equal to zero [3, [Il[ . Such an initial state has also 
been used in theoretical analysis of optomechanical sys- 
tems earlier [2^ . The state vector at a later time is then 
found to be 



|^(i))=^c(n,t)|^(n)). 



(19) 



n=0 



where \^{n)),n = S'[^(n)]|0)„i (with the 71-dependence of 
the squeezing parameter ^ emphasized) is the squeezed 
vacuum state of the mechanical oscillator, and 



c(n,t) 



^n -^[\a\^-iic^r,^t-v)] 



(20) 



D. Entanglement 

In this section we discuss bipartite entanglement be- 
tween the optical mode and the mechanical oscillator, 
which is of interest for quantum information processing 
as well as fundamental studies of quantum mechanics 

US US mi. 

A simple analytical understanding of entanglement in 
the present system may be arrived at by considering 
Eq. (|19p . The state \ip{t)) becomes disentangled when 
it can be written simply as a tensor product of the me- 
chanical and optical states. In order for this to occur 
either the optical or the mechanical state should factor 
out from the infinite sum of Eq. ([T^ . This can happen 
only if the mechanical state |^(n)) in each term is the 
same. This last condition may be satisfied if ^(n), the 
squeezing parameter, becomes independent of n. From 
Eqs. P^ - ([T5)) . it can be seen that the last factor in ^, 
i.e. sinh^ (?'sin;\;i/;\;), is independent of n only at i = 0, 
at which time ^(n) = 0. Thus, although the optical and 
mechanical modes are initially separable, they are entan- 
gled at all subsequent times. This behavior is in complete 
contrast to the case of linear optomechanical coupling, 
where in a similarly prepared system, the optical and 
mechanical modes disentangle once every mechanical pe- 
riod m. 

In our discussion of entanglement, we have neglected 
dissipation and noise. In the presence of these influences, 
the optomechanical states do not remain pure and be- 
come mixed, and entanglement is expected to weaken 
and eventually vanish as light ultimately escapes from 
the resonator. The entanglement of mixed states can be 
characterized by measures such as the logarithmic neg- 
ativity, which can be computed using, for example, the 
Quantum Langevin approach [32|. As this treatment is 
quite involved and not very revealing in the present case, 
we have not included it here. 



E. Phononic properties 

Knowledge of \ip{t)) from Eq. p^ allows us to inves- 
tigate the physical behavior of the mechanical oscillator, 
such as the average phonon number (b^b), whose time 
dependence is shown in the absence and presence of dis- 
sipation in Figs. HJa) and (b) respectively. 

We also show the temporal behavior of the phononic 
Mandel parameter [Eq. ([TU|] in Fig. [SJa). It can be 
seen from this figure that in the absence of dissipation, 
Qk = — 1 initially as expected for a vacuum state, and 
subsequently Q^ ^ 1 corresponding to a weakly squeezed 
vacuum state. The oscillations in Qk are due to the pe- 
riodic time dependence of {b^b) and {{b^b) ). The intro- 
duction of dissipation damps these oscillations as shown 




FIG. 4: Time evolution of tlie average phonon number (6^6) 
of the mechanical oscillator in the state Eq. (|19l) . (a) The 
phonon number evolution in the absence of noise and damping 
for g/ujm = 0.003, and a — 0.1. This curve has been scaled 
up by a factor of a thousand to improve visibility (b) The 
phonon number evolution obtained from the master equation 
of Eq. (|22p and including the effects of dissipation {'jo/oJm ~ 
0.1,7™/a;„ = lO^'^.n™ = 1000, T = 100 mK). It shows the 
absence of the coherent oscillations and an increase in the 
oscillator phonon number upon coming into contact with the 
thermal reservoir. At longer times, not shown here, {& &) 
reaches the steady state value {& b) — fim ~ 1000. 



in Fig. [5jb) and at longer times the value of Qk tends to 
the average number of bath phonons, as expected for a 
state in thermal equilibrium. 
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FIG. 5: The time dependence of the Mandel parameter [Qfc, 
defined in Eq. (|10p ] for the mechanical oscillator in the state 
Eq. p9[l . (a) The time evolution of the Mandel parameter for 
g/uim ~ 0.003, a = 0.1 in the absence of dissipation. From 
this curve it can be seen that at t = 0, the Mandel parame- 
ter is -1, corresponding to the occupation of the mechanical 
ground state. Later in time, Qh reaches the value 1, as ex- 
pected for a weakly squeezed vacuum state. The subsequent 
oscillations are due to the time variation in (6^b) and {{b^b)'^). 
(b) Time evolution of the Mandel parameter obtained from 
the master equation of Eq. (|22|) including the effects of dis- 
sipation {■yo/oJm = 0.1, ym/l^ni = 10^^,71™ = 1000, T = 100 
mK). The curve has been magnified by a factor of a hundred 
to provide better visibility, and shows the gradual disappear- 
ance of the coherent oscillations and the onset of increasing 
values oi Qk- At longer times, not shown in the plot, the 



Mandel parameter goes to the limit Qk 



as expected 



for a mechanical state in thermal equilibrium with the bath. 



F. Nonclassical state preparation 

In this section wc discuss the preparation of nonclas- 
sical states of mechanical motion, which arc of inter- 
est to metrology [20|, [23, [SJ]. We focus on mechanical 
states with sub-vacuum uncertainty (squeezing) in some 
quadrature. In general, this quadrature isa linear com- 
bination of position and momentum of the mechanical 
oscillator. We find that such squeezed states can be ob- 
tained from our system in two ways. 

The first technique can be understood from an inspec- 
tion of the reduced density matrix of the mechanical os- 
cillator. 



°'"E 
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2n 



-\an)){C{n)\, 



(21) 



which shows that the oscillator generally exists in a su- 
perposition of squeezed vacua. For a small number of 
average photons in the cavity (a^ ^ 1), only the first 
few terms in Ec^. (j21[) contribute to the quantum state. 
This can be verified by examining the time-evolution of 
the Wigner function corresponding to p,„, as shown in 



Fig. ini It can be seen that starting initially from the 
vacuum, the oscillator evolves generally into a state with 
sub-vacuum uncertainties in some quadrature. For the 
larg est available experimental coupling, {g/uJm ^ 0.003 
[12| ) we found, even in the absence of dissipation, very 
little squeezing. To obtain noticeable squeezing we in- 
creased the coupling to g/iOm = 0.08. In that case we 
found ~ 1.8 dB of squeezing at uJmt = 0.5, as shown 
in Fig. [6l In the presence of dissipation, as shown in 
Fig. [3 the squeezing degrades a little at uj,nt = 0.5, and 
more noticeably at w^i = 4. We note that unitary evo- 
lution in the quadratically coupled system yields super- 
positions of squeezed states rather than pure squeezed 
states. Nonetheless, this result is in complete contrast 
to the case of linear coupling, where unitary evolution 
cannot produce nonclassical mechanical states [20|. 

A more involved method can be used to produce pure 
squeezed states with quite large squeezing. In this ap- 
proach a larger number of average photons in the cavity 
(a^ ~ 25) is required, and in that case several terms in 
the sum of Eq. (j2ip make a non-negligible contribution to 
the reduced mechanical density matrix. The correspond- 
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FIG. 6: Time evolution of the Wigner function of the mechanical oscillator, found from the system wavevector in Eq. (|19p : q 
and p correspond to the position and momentum of the oscillator, respectively. The parameters used are g/iUm = 0.08, a = 1. 
No dissipation has been included. At ujmt — the initially occupied mechanical vacuum state can be seen. The squeezing at 
ijJmt — 0.5 is about 1.8 dB. 



ing Wigner function is plotted in Fig. |S1 and shows mul- 
tiple mechanical vacua differing in both the magnitude 
as well as the angle of their squeezing. The technique for 
producing pure squeezed states from this superposition 
follows from the form of the system state vector provided 
in Eq. (|19p , which shows that a conditional measurement 
of the photon number n in the cavity collapses the me- 
chanical state into the corresponding squeezed vacuum, 
\^{n))m [201 ■ This procedure can yield quite large me- 
chanical squeezing of about 10 dB. 



The large photon numbers in the cavity required for 
this technique, and the large phonon numbers resulting 
from mechanical squeezing unfortunately make the inclu- 
sion of dissipation effects computationally intensive. A 
treatment of the high photon number limit for our model 
as well as an extension to the regime of strong optome- 
chanical coupling will be the subject of a later work. We 
note that an analysis of mechanical squeezing at high 
photon numbers, using a different approach, has been 
presented in Ref . |15| . 



III. MASTER EQUATION TREATMENT OF 
DISSIPATION 

In this section we describe how we account for the pres- 
ence of environmental effects represented by Ho and Hm 
in Eq. ([1]). These terms can be handled in a prescriptive 
manner using the standard master equation approach 
[34| . The master equation is given by 

P = -t[p^Hq] + -^ {2apa'^ - a'^ap- pa^a) (22) 

+ ^ («>n + 1) (2fop6t - b^bp - pb^b) 

+ :(^n^ {2b^ pb - bb^ p - pbb^) , 



where 70 and 7^ are the decay rates of optical and me- 
chanical quanta to the respective reservoirs, and 



1 



,/fcBT _ 2' 



(23) 



the 



is the number of equilibrium quanta of energy hujr, 
mechanical bath which is held at temperature T. 

The master equation provided above is valid if g is 
smaller than either 70 or 7™. This condition ensures 
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FIG. 7: Time evolution of tlie Wigner function of the mechanical oscillator, obtained from the master equation [Eq. H22|) ]. The 
parameters used are g/ujm — 0.08, a = 1,'yo/ujm = 0.1,'ym/i-Om = 10~ ,nm = 1000, T = lOOmK. Comparing with Fig. |6]it can 
be seen that dissipation slightly changes the squeezing magnitude and angle at Umt = 0.5, and more noticeably at iUmt = 4. 



that the regime of strong optomechanical coupling which 
corresponds to g > (70, 7m), is avoided. This is be- 
cause in the strong coupling regime a master equation 
has to be derived in the dressed state basis in order 
to achieve consistency with thermodynamics |35| . We 
have also assumed that the electromagnetic vacuum is 
effectively at zero temperature. This is a reasonable ap- 
proximation since there are few thermal photons avail- 
able at optical wavelengths [13- Unless explicitly men- 
tioned, our calculations use experimentally accessible pa- 
rameters, g/ujm = 0.003, a -^ 1,70/0;™ = 0.1, 7™/^™ = 
10-7, n„ = 1000, and T = 100 mK [i,[i3. 

We have solved the master equation numerically in a 
frame rotating at the frequency wq for the initial state 
given by Eq. ^TE\\ . The transformation to the rotating 
frame removes uq from the dynamics, and thus we do 
not provide its value here. We have then obtained the 
system density matrix p as a function of time. All in- 
formation subsequently presented about the system has 
been obtained from the density matrix. The effects of 
dissipation have been discussed along with the analytical 
results for the dissipationless model in SecHIl We have 
verified that all our analytical results are reproduced by 
the master equation simulation in the absence of dissipa- 
tion. The effects of driving will be discussed in a future 



publication. 



IV. CONCLUSION 

We have explored the basic quantum characteristics 
of an important class of optomechanical system where 
the coupling varies as the square of the mechanical dis- 
placement. We have presented the dressed states, spec- 
trum, time evolution, entanglement, phonon statistics 
and wavefunction engineering relevant to the system. We 
have shown that in contrast to the case of linear cou- 
pling, quadratic interactions lead to more persistent en- 
tanglement, higher spectral nonlinearity, and a method 
of generating mechanical squeezing by unitary evolution. 
We have included dissipation in our analysis for the case 
of weak optomechanical coupling. Our results should 
be specifically of use to optomechanical spectroscopy, 
state transfer, wavefunction engineering, and entangle- 
ment generation. More generally, our work will be rele- 
vant to quantum information processing and sensing pos- 
sibilities opened up by optomechanical systems in the 
quantum regime. 
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FIG. 8: Time evolution of the Wigner function of the mechanical oscillator. The parameters used are g/uim ~ 0.08, a = 5. 
No dissipation has been included. At t = the initially occupied vacuum can be seen. At ujmt — 0.5, three ellipses, each 
corresponding to a squeezed mechanical vacuum appear. At cOmt = 9, four ellipses can be detected. 
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